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ABSTRACT 


A  new  composite  armor  concept  that  encompasses  an  extremely  hard  first  layer  to  deform 
the  projectile,  an  ortho  tropic  second  layer  to  slow  down  the  shock  wave  propagation,  a 
third  porous  layer  to  absorb  the  shock  wave  energy  through  PV-work,  and  a  fourth  layer 
to  provide  confinement  for  the  porous  medium  has  been  conceived.  Ceramic  Corbit-98 
and  Ceramet  Tungsten  Carbide  were  selected  as  candidate  materials  for  the  first  layer  and 
these  were  the  focus  of  this  research.  Dynamic  loading  responses  of  the  material  studied 
were  determined  through  planar  impact  experiment  conducted  on  a  single  stage  light-gas 
gun  at  NPS  Impact  Physic  Lab.  Impact  velocities  ranged  from  0.2  to  0.35  km/s.  The 
impactor  material  for  asymmetric  experiments  was  z-cut  single  crystal  sapphire. 
Diagnostics  used  included  a  VISAR  system,  to  measure  particle  velocities,  PZT  pins  to 
measure  onset  of  impact,  and  contact  pins  to  measure  impactor  velocities  and  tilt  angles. 
Through  this  study,  dynamic  loading  response  of  ceramic  Corbit-98  and  ceramet  tungsten 
carbide  were  determined.  The  Hugoniot  EOS  for  Corbit-98  ceramic  and  GC-915  tungsten 
carbide  were  measured  to  be  =  0.959  -i- 10.57  and  =  10.2  -i- 5.42  , 

respectively.  The  Hugoniot  Elastic  Eimit  (HEE)  of  GC-915  was  found  to  be  0.935  GPa 
and  spall  strength  of  approximately  2  GPa  was  also  measured. 
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I.  INTRODUCTION 


A,  BACKGROUND 

This  thesis  is  a  continuation  of  previous  Naval  Postgraduate  School  (NPS) 
research  performed  by  Denzel  (2010)  to  investigate  the  shock  properties  of  ceramic 
Corbit-98  (98%  Alumina)  and  above  that,  a  new  investigation  into  Tungsten-Carbide  in  a 
Cobalt  matrix  (WC/Co)  ceramet. 

The  main  motivation  of  this  research  is  based  upon  the  concept  of  employing  a 
layered  structure  (Ong,  2009)  for  personnel  armor.  This  layered  structure  stops  the 
projectile  in  a  series  of  stages  as  shown  graphically  in  Figure  1:  first,  a  projectile 
encounters  a  high  yield  strength,  high  impedance  material  which  causes  considerable 
plastic  deformation  of  the  projectile,  which  serves  twofold  purposes:  a)  to  transfer 
projectile  kinetic  energy  into  internal  energy  of  the  projectile  during  the  plastic  flow,  and 
b)  reduces  the  amplitude  of  the  shock  waves  produced  by  increasing  the  cross-sectional 
area  of  the  projectile;  second,  rapidly  spreading  the  shock  wave  created  by  the  projectile 
laterally  with  the  use  of  special  orthotropic  composite  material;  third,  converting  the 
remnants  of  the  kinetic  energy  from  the  impact  into  waste  heat  with  the  use  of  porous  or 
phase  changing  materials;  and  fourth,  a  final  stopping  layer  made  of  conventional  armor 
materials  such  as  Kevlar.  The  success  of  this  multi-layered  structure  for  personnel  armor 
will  allow  potential  space  and  weight  savings  due  to  the  associated  better  specific  (per 
unit  mass)  penetration  resistance  abilities.  The  focus  of  this  thesis,  is  on  the  first  layer 
(high  strength  and  high  impedance),  which  will  be  investigated  in  detail. 
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Figure  1.  Graphical  illustration  of  new  layered  armor  concept  (From  Ong,  2009) 

To  help  disperse  the  momentum  of  the  projectile,  it  needs  to  be  plastically 
deformed  causing  the  projectile  kinetic  energy  to  be  converted  into  internal  energy  and 
waste  heat.  In  order  for  that  to  happen,  the  armor  must  possess  a  much  higher  yield 
strength  so  thereby  causing  the  projectile  to  yield  before  the  armor.  This  yielding  point  is 
at  the  onset  of  plastic  deformation,  and  is  also  known  as  Hugoniot  Elastic  Limit  (FIEL)  in 
shock  physics. 

A  high  HEL  is  critical  for  any  armor  protection  material  to  cause  significant 
deformation  to  occur  in  the  projectile.  Therefore,  the  material  of  choice  for  the  first  layer 
would  be  a  material  of  very  high  HEL.  Due  to  the  high  HEL  of  Ceramic  Corbit-98  (8.27 
GPa),  it  was  selected  as  the  choice  of  material  to  be  studied.  Tungsten  Carbide  is  a  metal 
alloy  with  attractive  compressive  and  tensile  strength  properties,  making  it  ideally  suited 
for  use  as  a  protective  element  to  mitigate  shock-induced  effects  and  was  therefore 
chosen  as  the  second  material  to  study  for  the  purpose  of  this  thesis. 

The  prime  objective  of  this  research  is  to  experimentally  determine  the  HEL 
WC/Co  alloy  to  determine  their  suitability  as  the  Momentum  Dispersion  layer. 
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B. 


LITERATURE  RESEARCH 


1,  Ballistic  Protection 

At  present,  most  of  the  researeh  on  armor  is  based  on  finding  a  single  material 
that  ean  resist  penetration  from  projeetiles.  The  most  eommon  ehoiee  of  material  used  is 
Rolled  Homogeneous  Armor  (RHA),  whieh  is  a  type  of  hardened  steel  alloy  with  a  HEL 
of  between  2  to  3  GPa.  A  major  drawbaek  of  steel  is  its  high  density,  whieh  eauses  it  to 
be  an  impraetieal  eandidate  for  body  armor  due  to  the  weight.  However  with  new 
researeh  and  development,  the  material  eandidates  evolved  to  the  use  of  teohnieal 
eeramies  (e.g..  Aluminum  Oxide,  Boron  Carbide,  Aluminum  Nitrate  and  others  for  their 
high  strength  property  at  relatively  lightweight.  The  use  of  even  more  advanee  materials 
(e.g.,  Kevlar  Fiber-  Reinforeed  Polymers  [KFRP]),  Carbon  Fiber  Reinforeed  Polymers 
(CFRP),  and  Aramid  or  Polyethylene  woven  fabrie  eomposites  has  also  beeome  more 
eommon.  Figure  2  shows  a  brief  elassifieation  of  these  advaneed  eomposites. 
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Figure  2.  Relative  strength/tensile  modulus  of  advaneed  fiber  eomposites  (From  Ong, 

2009) 
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Such  evolution  of  protection  technology  has  had  varying  success  in  the  defeat  of  a  certain 
class  of  projectiles  and  it  became  evident  that  such  existing  technologies  may  have 
reached  a  plateau  in  the  development  of  body  armor. 

2.  Impetus  for  Ongoing  Research 

With  quickly  evolving  projectile  threats  and  stagnating  development  of  body 
armor,  the  need  for  new  armor  protection  concepts  to  break  through  the  current  paradigm 
becomes  clear.  Much  interest  has  been  shown  in  the  development  of  armor  protection 
using  layered  armor  as  demonstrated  by  Robbins,  Ding,  and  Gupta  (2004).  Gupta  and 
Ding  (2002)  demonstrated  the  use  of  wave  spreading  material  to  dissipate  the 
compressive  forces  of  the  incoming  projectile  within  a  microsecond  timescale.  Wilkins 
(1978)  have  also  shown  the  effectiveness  of  ceramics  in  plastically  deforming  the 
projectile  thus  defeating  it  from  the  onset  and  preventing  extensive  damage  to  the  lower 
layers  of  armor.  Herrman  (1969)  has  also  demonstrated  the  effectiveness  of  porous 
materials  in  absorbing  energy  during  shock  compression.  All  the  works  mentioned  above 
have  shown  promise  when  used  in  individual  layers,  but  there  is  still  a  need  to  put  these 
concepts  together  as  a  multi-layered  armor  system,  with  each  layer  playing  a  specific  role 
in  defeating  projectile  penetration. 

Poh  (2008)  has  shown  the  feasibility  of  composite  layered  armor  construction 
consisting  relatively  dissimilar  materials  with  each  material  aiding  to  resist  penetration 
with  specific  properties.  The  first  layer  of  the  composite  layered  armor  he  modeled 
consist  of  a  high  strength  layer  to  plastically  defeat  the  projectile,  and  a  wave  spreading 
idealized  material  as  a  second  layer  to  laterally  dissipate  the  compressive  shock  waves 
and  finally  a  porous  layer  to  assist  in  absorbing  the  residual  energy.  Through  numerical 
simulation  using  the  Autodyn®  hydrodynamic  code,  he  was  able  to  demonstrate  the 
feasibility  of  a  composite  layered  construction,  in  the  sequence  mentioned  above,  to  resist 
the  impact  of  a  projectile,  which  performed  better  than  an  AISI  4340  armor  grade  high 
strength  steel  plate  of  equivalent  thickness.  Figure  3  shows  results  from  his  simulations 
where  AISI  4340  steel  armor  is  perforated  as  compared  a  composite  armor  with  just 
minimal  penetration. 
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Figure  3.  Simulation  by  Poh  (2008)  to  demonstrate  capability  of  composite  armor  as 

compared  to  AISI  4340  steel. 


In  addition,  Ong  (2009)  has  conducted  live  firing  experiments  to  validate  this  new 
composite  armor  against  numerical  simulations.  Through  this  study,  the  composite  armor 
was  shown  both  experimentally  and  numerically  to  be  more  effective  in  resisting 
penetration  than  conventional  high  strength  armor  of  equivalent  (and  slightly  greater) 
areal  density,  and  that  the  material  layering  sequence  is  fundamentally  correct. 

In  order  to  reap  the  apparent  benefits  of  using  computers  to  simulate  the 
interactions  of  projectile  and  armor  systems,  the  description  of  materials  and  their 
response  within  the  extreme  terminal  ballistic  environment  must  be  accurate.  Although 
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the  simulation  mentioned  above  has  proven  the  effeetiveness  of  the  multi-layered  armor 
coneept,  sueh  simulation  have  not  been  able  to  preeisely  model  the  aetual  response  of  the 
material.  As  pointed  out  by  (Anderson  et  ah,  2009),  the  Drucker-Prager  method  fails  to 
aceurately  model  the  dynamic  response  of  ceramics  and  other  brittle  materials.  This 
further  reinforces  the  need  to  better  characterize  material  dynamic  response  in  order  to 
improve  the  accuracy  of  computational  simulations,  which  leads  to  the  main  objective  of 
this  thesis. 


3.  Research  Approach 

The  thesis  will  begin  with  an  overview  of  shock  physics  and  planar  impact 
experiments  before  going  into  the  experimental  setup.  As  this  is  a  continuation  of 
previous  research  of  Denzel  (2010),  the  research  will  continue  to  focus  on  validating  the 
dynamic  response  of  ceramics  as  well  as  a  new  investigation  into  ceramet  Tungsten 
Carbide.  Several  single  stage  light-gas  gun  experiments  have  been  performed  to  acquire 
dynamic  data  that  can  better  characterize  the  material  model  used  for  computational 
simulations.  Finally  experimental  results,  data  analysis  and  final  technical  results  will  be 
discussed  and  recommendation  will  be  made  for  future  research. 
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II.  SHOCKWAVE  THEORY 


A  good  understanding  of  the  dynamic  response  of  materials  studied  under  rapid 
impulsive  loading  is  vital  to  the  use  of  these  materials  in  applications  involving  high 
velocity  impact.  Therefore  there  is  an  impetus  to  identify  the  strength  properties  and  the 
inelastic  deformation  mechanisms  of  the  chosen  material  under  well-characterized,  planar 
shockwave  compression.  This  section  describes  the  literature  research  done  on  the 
fundamentals  of  shock  physics  that  will  be  applied  to  the  experiments  conducted  for  the 
purpose  of  this  thesis. 

A,  SHOCKWAVE  FUNDAMENTALS 


Strain  produced  in  materials  is  directly  proportional  to  the  stress  applied  to  it. 
This  linear  behavior  only  holds  true  up  to  a  point  at  which  the  material  will  not  return  to 
its  exact  original  shape  when  stress  is  released.  This  point  is  known  as  the  yield  point  or 
elastic  limit  of  the  material.  When  a  material  is  strained  beyond  its  elastic  limit,  plastic 
deformation  occurs.  We  can  see  how  this  affects  material  response  to  shock  compression 
to  high  pressure  by  looking  at  some  fundamental  concepts. 


Sound  velocity  is  constant  in  the  elastic  region  of  the  material  and  guided  by  the 
following  equation: 


dp  ' 

This  implies  that  in  the  elastic  region,  pressure  and  density 
the  elastic  region,  P/  p  is  not  linearly  proportional  and 
pressure  or  density.  This  nonlinearity  can  be  observed  from 


(1) 

are  linearly  related.  Beyond 
sound  speed  increases  with 
Figure  4. 
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A 


Figure  4.  Amplitude  wave  in  regime  of  pressure  beyond  yield  strength  (From 

Wikibooks,  2009) 

Due  to  the  nonlinearity  observed,  point  C  of  the  waveform  will  have  the  lowest  local 
wave  speed  while  points  B  and  A  will  have  a  higher  wave  velocity  due  to  higher 
pressure.  This  difference  in  velocity  will  eventually  allow  point  B  and  A  to  catch  up  as 
depicted  by  the  time  lapse  diagram  illustrated  in  Figure  5. 


Figure  5.  Steepening  of  pressure  wave  (From  Wikibooks,  2009) 
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With  time,  the  pressure  wave  will  eventually  steepen  at  the  front  to  form  a  vertieal  line 
that  is  known  as  a  shoekwave.  This  oeeurs  beeause  for  normal  materials,  sound  speed 
inereases  with  inereasing  pressure. 

B,  DISCONTINUITY 

When  a  pressure  wave  steepens  to  beeome  a  shoekwave,  a  large  gradient  in 
temperature  and  veloeity  ean  be  observed  at  the  front  of  the  shoek.  This  abrupt  ehange 
from  the  non-shoeked  to  the  shoeked  stated  is  often  taken  to  be  a  diseontinuity. 

On  a  mieroseopic  level,  large  gradient  in  temperature  and  veloeity  at  the  front  of 
the  steepening  wave  cause  an  irreversible  process  to  occur  which  cause  entropy  to 
increase.  The  net  effect  on  a  macroscopic  level  is  that  mass,  momentum,  and  energy  is 
conserved  across  the  shock  front,  but  entropy  is  not.  Rankine,  Rayleigh,  and  Hugoniot 
showed  that  an  adiabatic  shock  front  would  violate  conservation  of  energy,  and  therefore 
shock  fronts  must  be  non-adiabatic  and  irreversible  (Los  Alamos  Science,  1985). 


Figure  6.  Irreversible  changes  by  shock  front. 

From  Figure  6,  we  can  see  that  the  internal  energy  increase  is  proportional  to  the 
area  under  the  Rayleigh  line.  If  we  take  the  release  path  to  be  close  to  the  Hugoniot,  the 
material  will  follow  this  curve  down  when  pressure  is  released  and  result  in  some  residual 
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energy  (heat)  left  in  the  material.  Therefore,  it  is  important  to  note  that  the  shoek  proeess 
is  irreversible  and  the  material  is  ehanged  as  a  result  of  having  been  shoek  eompressed. 
For  shook  loading  to  a  very  high  pressure,  the  residual  heat  oan  be  substantial. 

C.  RANKINE-HUGONIOT  JUMP  CONDITIONS 

Due  to  reasons  disoussed  above,  the  original  states  of  partiole  velooity  Uo,  density 
po,  internal  energy  Eq  and  pressure  Pq,  suddenly  ohanges  aoross  the  shook  front.  As 
mentioned  above,  mass,  momentum,  and  energy  are  oonserved  aoross  the  shook  front. 
The  oonservation  relations  do  not  depend  upon  a  proeess  path  but  merely  upon  the  initial 
and  final  states  of  the  material  in  question.  Three  equation  are  oommonly  derived  from 
the  oonservation  rules,  and  are  known  as  “Rankine-Hugoniot  Jump  equations”.  Based  on 
a  Lagrangian  ooordinate  system,  the  jump  equations  are  expressed  as  shown  below: 

Mass:  ^  =  1  -  ^  (2) 

P  Us 

Momentum:  P-  P^=  pJJJU ^  (3) 

Energy:  E  -E^  =  —(P  +  PJ(V^  -  V)  Where  V=l//?  (4) 

With  3  equation  and  5  unknowns  (e.g..  Us,  Up,  P,  and  E),  it  beoomes  olear  that 
another  relationship  is  required  to  relate  some  of  the  states.  This  equation  is  often 
referred  to  as  Equation  of  state  (EOS).  Unfortunately,  there  exists  no  EOS’s  that  oan  be 
derived  from  first  prinoiples  to  desoribe  non-linear  elastio  materials.  This  drives  us  to 
measure  two  of  the  five  unknown  quantities  experimentally.  This  in  turn  allows  us  to 
define  an  EOS. 

D.  SHOCK  HUGONIOT  IN  Us  -  Up  PLANE 

In  the  past,  many  experiments  were  oonduoted  to  determine  EOS’s  for  many 
materials,  and  it  was  found  that  shook  and  partiole  velooity  are  linearly  related  in  most 
materials  by  the  equation  shown  below: 

U,  =  c„  +  su^ 
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Therefore  if  Us  and  Up  are  determined  experimentally  at  a  number  of  different 
shoek  states  for  a  certain  material,  the  shock  Hugoniot  in  Us  -  Up  plane  can  then  be  easily 
described.  When  the  data  points  are  plotted  in  the  Us  -  Up  plane,  doing  a  linear  fit 
through  the  points  would  give  the  constant  Co  at  the  y-axis  intercept,  and  also  determine 
the  slope.  The  bulk  sound  speed  Cb  is  typically  a  very  good  approximation  for  this 
constant  Cq.  Figure  7  illustrates  an  example  of  a  Us  -  Up  Flugoniot  plot  for  Tungsten 
Carbide. 


Figure  7.  Hugoniot  lot  in  Us  -  U?  of  tungsten  carbide  (From  Grady,  1995) 

Hence,  an  EOS-derived  Hugoniot  can  then  be  determined  from  experimental  data.  The 
acquisition  of  these  experimental  data  is  exactly  the  main  focus  of  this  thesis  research  in 
order  to  fully  describe  the  shock  states. 

E.  PLANAR  IMPACT  EXPERIMENTS 

Light-gas  guns  are  commonly  used  to  launch  flat-plate  impactors  onto  flat  plate 
targets  to  produce  normal  shocks  into  the  targets.  These  normal  shocks  are  critical  to 
produce  the  uni-axial  strain  conditions  required  to  use  the  simple  theories  described 
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above  and  thus  obtain  EOS  information  from  simple  particle  velocity  and  shock  speed 
measurements.  This  method  offers  tight  control  over  the  stress-time  history  of  the 
sample.  Principle  diagnostics  are  the  Laser  Velocity  Interferometer  for  Any  Reflector 
(VISAR),  and  time  of  arrival  pins. 

The  experiment  usually  involves  a  projectile  impacting  a  stationary  target  plate  at 
the  velocity  of  Ud.  A  discontinuity  of  particle  velocity  and  stress  between  the  target  and 
projectile  will  be  formed  on  impact,  and  shockwaves  are  immediately  developed  to 
relieve  this  discontinuity.  The  developed  shockwaves  will  then  propagate  away  from  the 
plane  of  impact  in  both  the  projectile  and  target.  With  the  help  of  the  VISAR  system  and 
PZT  pins,  the  velocity  of  the  shockwave  can  then  be  determined.  Projectile  velocities  Ud 
is  measured  by  shorting  6  electrically-charged  pins  located  at  measured  distances 
(millimeters  ahead  of  the  target  sample)  and  subsequently  particle  velocity  Up  can  be 
either  measured  with  the  VISAR  system  or  determined  using  the  jump  conditions.  These 
parameters,  together  with  the  help  of  fundamental  equations  will  allow  us  to  develop  an 
EOS  under  dynamic  loading  conditions. 

As  mentioned  above,  to  determine  the  fundamental  shock  response  of  a  target 
material,  a  plot  of  shock  wave  velocity  Us  versus  particle  velocity  Up  is  required. 
Therefore  several  impact  experiments  at  various  Ud  will  have  to  be  conducted. 
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III.  EXPERIMENTAL  METHODS 


The  measurement  of  dynamic  properties  of  materials  requires  the  use  of  well- 
developed  experimental  techniques  and  often  sophisticated  equipment  to  obtain 
satisfactory  data.  In  this  section  of  the  thesis,  equipment  used  and  theory  behind  how  the 
experiment  is  conducted  to  acquire  the  data  will  be  elaborated. 

A.  LIGHT-GAS  GUN 


Figure  8.  NFS  light-gas  gun  assembly 


All  the  experiments  conducted  for  this  thesis  was  done  on  a  76  mm  bore  single 
stage  light  gas  gun  housed  in  the  impact  physic  lab  of  the  Naval  Postgraduate  school 
(NFS)  as  shown  in  Figure  8.  This  helium-driven  gas  gun  was  designed  to  accelerate 
projectiles  up  to  velocities  below  0.5  km/s,  thus  creating  high-pressure  states  that  can  be 
used  to  characterize  the  fundamental  dynamic  responses  of  the  material  tested.  The 
breech  utilizes  a  “wrap  around”  design  that  eliminates  the  need  of  any  diaphragms. 
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High-pressure  helium  gas  eontained  in  the  breeeh  is  separated  from  the  launeh  tube  with 
the  use  of  o-rings  installed  on  both  ends  of  the  projeetile.  Figure  9  shows  the  pieture  of 
the  projeetile  and  the  breeeh: 


Breeeh 

Assembly 


Inlet 

(Vaeuum/Low 

pressure) 


High  Pressure 
Chamber 


Launeh  Tube 


Figure  9.  Breeeh  assembly  and  projeetile  of  gas  gun  (From  Ho,  2009) 


The  projeetile  is  launehed  when  a  small  amount  of  low-pressure  is  introdueed 
behind  the  projeetile,  pushing  it  past  the  ports  of  the  breeeh  and  thus  allowing  the  high- 
pressure  helium  gas  to  aeeelerate  the  projeetile. 

To  ensure  that  the  launeh  tube  is  properly  aligned  to  the  target,  adjustments  were 
regularly  made  to  the  launeh  tube  support  assembly  with  the  help  of  a  elass  3b  He:Ne 
laser  sent  down  the  entire  length  of  the  launeh  tube  and  retro-refleeted  baek.  In  order  to 
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ensure  a  very  flat  impact  condition,  the  alignment  of  the  launch  tube  must  be  less  than  1 
mrad  relative  to  the  target  plane.  The  alignment  of  the  barrel  is  adjusted  as  shown  in 
Figure  10. 


Launch  Tube 
Support 


Figure  10.  Figure  of  launch  tube  support  assembly  (From  Ho,  2009) 

Finally  the  projectile  will  impact  onto  the  target,  that  is  mounted  in  diagnostic 
tank  assembly.  The  effect  of  air  buildup  between  the  target  and  projectile  is  minimized 
by  creating  a  vacuum  condition  of  <30  millitorr  in  the  barrel  and  the  diagnostic  tank 
assembly  prior  to  firing.  An  expansion  chamber  between  the  end  of  the  launch  tube  and 
target  assembly  is  provided  to  further  ensure  that  the  air  buildup  is  kept  to  a  minimum. 

Lastly,  the  catcher  tank  assembly  ensures  that  after  impact,  momentum  of  the 
debris  is  stopped  to  minimize  damage  to  the  gun.  This  is  achieved  with  several  layers  for 
energy  absorbing  material  in  the  form  of  aluminum  honeycomb  sheets,  a  blast  shield — an 
extra  layer  of  mild  steel  to  prevent  any  debris  from  striking  the  rear  surface  of  the  catch 
tank  and  a  sliding  baffle  to  absorb  energy.  The  baffle  has  nylon  wheels  and  is  a  sealed 
cylinder  of  steel.  To  ensure  that  the  honeycomb  can  be  uniformly  compressed,  a  thick 
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aluminum  place  is  placed  in  front  of  the  honeycomb  so  that  on  impact,  the  plate 
compresses  the  honeycomb  material  relatively  uniformly  causing  energy  to  be  dissipated. 
Figure  1 1  shows  the  catcher  tank  assembly. 
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Figure  1 1 .  Catcher  tank  assembly  (From  Ho,  2009) 


B.  SAMPLE  CHARACTERIZATION 

In  order  to  make  aceurate  dynamic  property  measurements,  we  must  first 
characterize  certain  initial  properties  for  the  selected  material  to  be  experimented.  This 
initial  characterization  is  important  to  establish  the  baseline  properties  for  the  material 
being  studied  and  later  to  be  used  for  the  aecurate  derivation  of  dynamie  loading 
properties  of  the  material.  Espeeially  for  sintered  materials  like  polycrystalline  ceramics, 
it  must  be  ensured  that  it  is  close  to  being  isotropic.  As  shown  by  Gust  and  Royce 
(1971),  previous  experimental  results  for  similar  materials  have  elearly  shown  that  initial 
sound  speed  depends  strongly  upon  initial  density  and  that  even  small  deviations  in  initial 
density  cause  measurable  differenee  in  initial  sound  speed  and  dynamic  properties. 
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Figure  12.  Longitudinal  sound  speed  vs.  density  of  eeramies  (From  Denzel,  2010). 

From  the  plot  in  Figure  12,  it  ean  be  seen  that  the  measured  density  and 
longitudinal  sound  speed  of  the  Corbit-98  is  in  agreement  with  the  previous  work  of  Gust 
and  Royce  (1971). 

1,  Initial  Density 

As  mentioned  in  Denzel  (2010),  the  initial  density  of  the  Corbit-98  sample  used  in 
the  experiments  of  this  thesis  was  measured  at  Los  Alamos  National  Laboratory  (LANL) 
using  standard  immersion  teehnique  that  is  based  on  Arehimedes  prineiple.  The  sample’s 
weight  is  measured  in  dry  air  and  submerged  in  water.  By  knowing  the  exact  temperature 
of  the  distilled  water,  one  can  obtain  the  volume  displaced  by  the  sample  by  dividing  the 
difference  in  weight  between  the  dry  and  the  underwater  measurements  by  the  water 
density.  Then,  knowing  the  samples  mass  and  volume,  once  can  very  accurately 
calculate  its  density.  The  density  was  found  to  be  =  3.864  +  0.005  g/cc. 

Density  of  the  Ceramat  Tungsten  Carbide  samples  was  determined  by  using  a 
mircometer  to  measure  diameter  and  thickness  and  an  electronic  scale  with  an  accuracy 
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of  ±0.1  gram  to  measure  their  weight.  This  is  possible  due  to  the  high  tolerance  for  the 
samples  as  supplied  by  General  Carbide,  the  density  was  found  to  be  14.03  g/cm  for 
samples  marked  GC-915  and  12.82  g/cm  for  samples  marked  GC-330,  which  coincide 
with  the  figures  published  on  the  General  Carbide’s  specification  sheet  for  the  material. 


2.  Elastic  Sound  Speeds 


Using  Hooke’s  law  along  with  Newton’s  Second  Law,  one  can  show  that  the 
speed  of  sound  within  a  material  is  a  function  of  the  innate  properties  of  the  material. 
The  general  relationship  between  the  sound  speed  in  a  solid,  density  and  elastic  constants 
is  shown  by  the  following  general  equation; 


C, 


(6) 


where  Ci  is  the  sound  speed,  M  is  the  appropriate  elastic  constant  and  p  is  the  material 
density.  It  must  be  noted  that  the  subscript  ij  of  the  elastic  constant  M  is  used  to  indicate 
the  wave  type  being  considered.  Types  of  possible  sound  wave  include  longitudinal 
sound  speed  (Cl),  shear  sound  speed  (Cs),  bulk  sound  speed  (Cb): 

Bulk  Sound  Speed:  -^C^;  (7) 


From  these  sound  speeds,  modulii  and  Poisson  ration  can  then  be  determined: 


Longitudinal  Modulus: 

II 

(8) 

Shear  Modulus: 

G  =  PoCl\ 

(9) 

Bulk  Modulus: 

II 

(10) 

i3K-F) 

Poisson  Ratio:  v  = - 

{F-2G) 

(3K-2G) 

(11) 

(3^  +  F) 

~  2(F-G)~ 

2(3K  +  Gy 

Both  Cl  and  Cs  of  the  sample  used  for  the  purpose  of  this  thesis  were  measured 
using  a  commercially  available  ultrasonic  pulse-echo  system  made  by  Panametrics  Inc.  (a 
division  of  Olympus).  This  system  consists  of  a  pulse  unit  that  can  be  used  with  either 
longitudinal  or  shear  transducers,  in  pulse-echo  geometry.  The  observation  and 
measurement  of  timing  between  pulses  were  done  using  a  fast  digital  oscilloscope.  More 
than  1  sample  of  each  material  was  measured  to  ensure  that  the  samples  were  of  uniform 
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quality.  For  each  sample,  many  reflections  were  measured  and  a  least  square’s  fit  was 
used  to  obtain  the  sound  speeds.  Table  1  summarizes  the  acoustic  properties  measured 
and  calculated  for  both  ceramic  Corbit-98  and  ceramet  tungsten  carbide; 


Table  1.  Properties  of  target  sample 


Properties 

Corbit  98 

WC/Co  GC-330 

WC/Co  GC-915 

Density  (q/cm^) 

3.864 

12.82 

14.03 

Elastic  Wave  Velocities  (km/s) 

Longitudinal 

10.55 

6.29 

6.67 

Shear 

6.18 

3.63 

3.92 

Bulk 

7.77 

4.69 

4.90 

Shear  Modulus  (GPa) 

148 

169 

216 

Longitudinal  Modulus  (GPa) 

430 

507 

624 

Bulk  Modulus  (GPa) 

233 

282 

337 

Poisson  Ratio 

0.239 

0.250 

0.236 

C.  HUGONIOT  MEASUREMENT  EXPERIMENTS 


As  mentioned  earlier,  one  of  the  main  objectives  of  this  thesis  is  to  perform 
impact  experiments  to  measure  shock  and  particle  velocity  over  a  range  of  impact  in 
order  to  determine  the  Hugoniot-EOS  for  the  material  studied.  In  this  section,  the  theory 
and  physical  set  up  of  the  gas  gun  planar  impact  experiment  will  be  explained  in  detail. 
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1,  Flyer  Velocity 


Figure  13.  Gas  gun  performance  plot 


By  design,  the  gas  gun  achieves  desired  projectile  speed  by  adjusting  the  breech 
pressure  based  on  the  projectile  weight.  This  calculation  is  based  on  the  gun  performance 
curve  shown  in  Figure  13.  After  the  projectile  is  fired,  its  velocity  before  impact  is 
accurately  determined  by  using  a  stepped  circular  array  of  electrical  shorting  pins 
mounted  on  the  target  plate.  The  main  purpose  for  accurately  measuring  this  velocity  is 
to  derive  the  particle  velocity  in  the  target,  which  will  be  elaborated  on  later  in  this 
section. 


2,  Shock  Speed  Us 

For  material  under  elastic  compression,  the  wave  speed  must  be  close  to  the 
longitudinal  sound  speed.  However  for  finite  amplitude  waves  where  shock-up  is 
expected,  shock  velocities  will  be  greater  than  Cb.  For  these  kinds  of  planar  experiments. 
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shock  speed  Us  is  aetually  not  direetly  measured.  It  is  in  faet,  derived  by  dividing  the 
thiekness  of  the  target  and  the  shoek  transit  time  as  shown  by  the  equation; 


U.  = 


X 

t 


Shook  transit  time  is  known  by  measuring  the  time  from  impaot  until  the  breakout  of  the 
first  wave  at  the  baok  of  the  target.  To  be  able  to  measure  this  short  transit  time  (on  the 
order  of  a  mioroseoond),  time  must  be  resolved  down  to  nanoseoonds.  This  is  aohieved 
through  the  use  of  piezoeleotrio  (PZT)  pins  and  interferometers. 


1  Trigger 


Velocity 

Interferometer 

(VISAR) 


! 


Velocity  Pin 


PZT 

Pins 


in 


Figure  14.  Sohematio  of  planar  impaot  setup 

From  the  diagram  shown  in  Figure  14,  it  oan  be  seen  that  a  trigger  pin  is  affixed  to 
the  target  so  as  to  trigger  the  eleotronio  system  upon  impaot  of  the  projeotile.  A  Velooity 
Interferometer  System  for  Any  Refleotor  (VISAR)  on  the  other  end  of  the  target  is  used 
to  measure  the  time  at  whioh  the  first  breakout  of  the  shookwave  appears.  Additional 
PZT  pins  allow  projeotile  impaot  time  to  be  measured.  As  some  materials  are  not 

21 


reflective  enough  to  allow  direct  interface  with  the  VISAR,  a  very  thin  foil  may  be 
applied  with  epoxy  to  the  target  to  increase  reflectivity.  Finally,  impact  velocities  were 
measured  by  shorting  electrically-charge  velocity  pins  located  at  measured  distances  a 
few  millimeters  ahead  of  the  target  sample. 

3.  Particle  Velocity 

For  planar  impact  experiments,  particle  velocity  is  closely  related  to  flyer 
velocity.  For  symmetrical  impact,  where  both  the  flyer  and  target  are  of  the  same 
material,  the  math  is  greatly  simplified  and  the  particle  velocity  can  be  taken  as  0.5 *Ud 
or  half  the  flyer  velocity. 

For  the  unsymmetrical  impact  experiment,  to  get  Up2  or  particle  velocity  in  the 
target,  we  use  the  fact  that  across  the  impact  interface,  pressure  and  particle  velocity  must 
be  conserved  as  shown  in  Figure  15. 


Figure  15.  x-t  diagram  for  unsymmetrical  impact 


Usually  for  the  an  unsymmetrical  impact  experiment,  Hugoniot  of  flyer  is  known 
but  flyer  velocity  and  shock  velocity  in  the  target  are  required  to  be  measured. 


For  the  flyer  (material  1),  in  general  we  have: 

Since  =U^-  ,  we  can  rewrite  Equation  (13)  as: 


(12) 

(13) 
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P,=p,, {S, (U^  -U,,)  +  Q  -  U,, ) 

=>  Pj  =  {Pq^S^)Ujj  —  2Pq^S^UjjU p2  +  Poi^l^Pl  +  Poi^J-^D  ~  Poi^J-^P2 

^P,=  (PoAWl^  -  (2AAt/z,  +  pMU,2  +  P,^(S,Ul  +  C,U^)  (14) 

Since  Us  can  be  derived  from  the  experiment,  pressure  P2  can  then  be  determined  by  the 
following  equation: 

P2=PoUs2Up2  (15) 

Sinee  pressure  Pi  and  P2  must  be  equal  (from  eonservation  of  mass  and  momentum),  we 
obtain 

(Poi^l)^ p2  ~  (^Poi^J-^D  PoiP'l)^P2  Poi(‘5'i1^£I  +  ~  Po2^S2^ p2 

(PoA)f^p2  -{^PoX^Pd  +P01Q  +  Pq2Us2)Up2  +  Po\{Sp\  +  Cjf/^,)  =  0  (16) 

From  Equation  (16),  Up2  can  then  be  found. 

D,  HUGONIOT  ELASTIC  LIMIT 

In  addition  to  establishing  a  point  on  the  shoek  Hugoniot  Us  -  Up  plane,  an 
important  objeetive  of  this  thesis  is  to  determine  the  dynamic  yield  point  or  the  Hugoniot 
Elastie  Eimit  (HEE)  qhel  for  the  materials  studied. 
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Figure  16.  Shock  Hugoniot  in  P-V  plane  for  a  material  with  strength 


From  Figure  16,  the  HEL  point  is  shown  at  the  end  of  the  slope  of  the  elastic 
wave,  as  an  inflection  point.  Below  the  HEL  would  be  the  elastic  region  and  above,  the 
plastic  region.  In  the  elastic  region,  the  velocity  of  the  wave  would  be  the  longitudinal 
sound  speed  Cl.  In  the  plastic  region,  the  slope  of  the  Raleigh  line  would  be  less  than 
that  of  the  elastic  Hugoniont  for  low  stresses,  and  this  results  in  a  shock  velocity  below 
Cl.  Eor  high  stresses  the  shock  speed  can  be  greater  than  the  longitudinal  sound  speed. 
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Figure  17.  A  typical  wave  profile  of  a  free-surface  VISAR  measurement  of  aluminum 

(From  Isbell,  2005) 


Due  to  the  different  sound  speeds  in  these  regions,  a  cusp  can  easily  be  observed 
in  the  wave  profile  when  the  material  is  stricken  up  to  the  plastic  region  as  illustrated  in 
Figure  17.  Therefore,  in  order  to  determine  Qhel,  first  of  all  the  amplitude  of  the  wave  at 
the  cusp  is  measured  while  taking  note  that  it’s  the  free  surface  velocity.  Since  the  wave 
would  be  travelling  at  the  longitudinal  sound  speed  Cl  at  the  yield  point,  and  the  initial 
density  of  the  material  can  be  measured,  qhel  can  then  be  derived  by  applying  the 
momentum  jump  condition: 

=  (17) 

For  an  elastic  wave.  Us  is  simply  Cl  and  hence, 

(^HEL=\iPoCLUl)  (18) 

The  above  equation  is  only  applicable  to  a  free  surface  experiment  where  no  window  was 
used. 
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E. 


SPALL  MEASUREMENTS 


Another  objective  of  this  thesis  is  to  measure  the  spall  strength  of  the  materials 
studied.  Spall  fracture  is  usually  a  result  of  the  interaction  of  compressive  stress  wave 
from  relatively  low-impedance  interfaces  (normally  free  surfaces)  and  subsequent  wave 
interactions  as  shown  in  Figure  18.  The  release  wave  interactions  cause  tension  to  be 
generated  in  the  target. 


Figure  18.  Interaction  of  release  wave 

The  spall  fracture  data  in  this  thesis  were  obtained  by  planar  impact  experiments 
where  uniaxial  strain  conditions  exist.  Such  an  impact  experiment  is  usually  designed  to 
allow  two  release  waves  interact  in  the  target  and  thus  bringing  the  material  into  tension 
dynamically.  If  the  tension  exceeds  the  spall  strength  of  the  target,  the  target  will  fail  and 
part  at  this  plane.  In  the  design  of  the  experiment,  the  thickness  of  the  target  and 
impactor  are  chosen  in  such  a  way  that  allows  the  release  wave  from  both  free  surfaces  to 
meet  at  approximately  at  the  center  of  the  target. 
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Figure  19.  A  typical  free  surface  VISAR  wave  profile  showing  a  spall  signature  (From 

Denzel,  2010) 

Figure  19  shows  a  typical  wave  profile  of  a  free  surface  planar  impact  experiment 
exhibiting  a  spall  signature.  The  simplest  but  also  a  least  accurate  way  to  estimate  the 
spall  strength  is  to  find  the  difference  in  peak  particle  velocity  and  minimum  particle 
velocity.  Then  applying  the  momentum  jump  condition  yields  the  following  equation; 

=  (19) 


T - 1 - 1 - 1 - 1 - r 
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The  wave  speed  Ux  used  in  a  duetile  material  shoeked  above  HEL  is  the  bulk  sound  speed 
(Cb).  As  for  brittle  materials  like  eeramies,  spall  strength  eannot  be  aeeurately 
determined  with  the  use  of  bulk  sound  speed  as  the  wave  speed  in  the  above  equation. 
However,  for  simplifieation  and  for  just  an  estimation  of  spall  strength,  bulk  sound  speed 
is  knowingly  used  to  aequire  an  approximate  value  of  the  spall  strength. 

F.  TARGET  SETUP 

1.  Target  Holding  Plate 

The  most  generie  target  holding  plate  used  for  the  purpose  of  this  researeh  is  a 
six-ineh  diameter  aluminum  6061  dise  as  shown  in  Figure  20: 


Figure  20.  Sehematie  of  target  holding  plate 

VISAR,  PZT  and  veloeity  pin  diagnosties  were  utilized  for  the  purpose  of  this 
researeh.  From  Figure  20,  it  ean  be  noted  that  the  target  holding  plate  has  a  total  of  7 
holes  for  veloeity  pins  (1  being  the  grounding  pin),  3  holes  for  PZT  pins,  one  for  the 

mounting  of  the  VISAR  probe  and  3  for  mounting  of  the  plate  on  to  the  diagnostie  tank. 
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Before  mounting  any  of  the  diagnostics  and  target  material,  the  plate  is  first  lapped  on  a 
fine  sand  paper  on  a  granite  flat  to  a  surface  finish  of  12  microns.  Similarly,  the  target 
sample  must  be  lapped  to  be  flat  and  parrallel. 

2,  Projectile 


Figure  21.  Aluminum  alloy  6061  projectile 

The  generic  projectile  as  shown  in  Figure  21  is  used  for  all  the  experiments 
conducted.  The  projectile  is  made  of  aluminum  alloy  6061  and  machined  to  a  tolerance 
of  less  than  one  thousands  of  an  inch  on  the  outer  diameter  as  well  as  on  the  impacting 
face. 


3.  VISAR 

The  laser  velocity  interferometer  measures  surface  or  interface  velocity  as  a 
function  of  time.  Either  free  surface  motion  or  motion  of  the  interface  between  the  target 
sample  and  a  transparent  window  material  can  be  determined  with  great  time  and  velocity 
resolution.  The  type  of  interferometer  used  for  the  purpose  of  this  thesis  is  a  velocity 
interferometer  for  any  reflector.  VISAR  is  a  modified  Michelson  interferometer  where 
velocity  instead  of  displacement  are  directly  obtained.  This  state  of  the  art  VISAR 
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system  was  acquired  from  National  Security  Technologies  (NSTEC).  With  a  small 
physical  footprint  and  ease  of  use,  this  system  makes  it  perfect  for  research  conducted  by 
students.  The  laser  utilized  for  this  system  is  a  class  4  frequency  doubled  NdiYAG  solid 
state  continuous  wave  (CW)  laser.  Procured  from  Coherent,  this  laser  is  designed  to 
operate  in  the  wavelength  of  532  nm.  The  delivery  system  for  the  laser  to  the  diagnostic 
tank  is  an  optical-fiber  system  that  is  safe  and  effective  method  for  coupling  the  laser 
light  into  the  target  chamber.  The  output  of  the  laser  is  focused  into  this  50  pm  step 
index  optical  fiber  through  voltage  controlled  wave  plates  (Pockels  cell).  An  optical  fiber 
probe  is  connected  at  the  end  of  this  optical-fiber  system  through  a  vacuum  flange  in  the 
diagnostic  tank.  The  laser  light  is  projected  onto  the  target  through  a  Plano-convex  lens 
affixed  at  the  tip  of  the  probe,  and  this  same  lens  serves  to  collect  the  reflected  light.  The 
reflected  laser  light  is  then  sent  back  to  the  interferometer  cavity  and  interference  occurs. 
Quadrature  signals  that  contain  information  about  the  interference  pattern  are  routed  to 
four  928  photomultipliers  before  the  converted  electrical  signals  are  sent  to  a  Tektronix 
DPO4104  transient  recorder.  The  recorded  data  is  then  finally  stored  on  a  connected 
laptop  ready  for  analysis. 

4.  Velocity  Pins 

As  mentioned  earlier  in  this  section,  there’s  a  need  to  measure  the  projectile 
velocity  and  this  is  achieved  with  a  stepped  circular  array  of  electrical  shorting  pins  as 
shown  in  Figure  22.  The  use  of  these  pins  allows  the  speed  as  well  as  the  angle  of  tilt  of 
the  projectile  to  be  measured.  A  total  of  seven  pins  are  used  for  every  experiment 
conducted  for  the  purpose  of  this  thesis.  Six  of  the  pins  are  aligned  at  60  degrees  apart 
from  each  other  in  a  circular  array.  The  7th  pin  or  the  ground  pin  is  placed  in  such  a  way 
that  it  will  be  the  first  to  be  impacted.  Every  pin  is  individually  shrink-wrapped  so  as  to 
be  electrically  isolated  from  the  target  holding  plate  till  the  impact  of  the  projectile 
occurs.  The  non-impact  side  of  the  pins  are  soldered  to  a  wirings  harness  that  is 
connected  to  a  pin  circuit.  This  pin  circuit  generates  an  electrical  pulse  for  each  pin  that 
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is  shorted  on  impact  by  the  projectile.  Using  a  least-squares  fit  to  the  times  at  which  the 
pins  were  shorted,  combined  with  the  various  lengths  and  angles,  both  projectile  velocity 
and  tilt  can  then  be  derived. 


5,  Piezoelectric  (PZT)  Pins 

In  order  to  determine  the  time  of  impact,  PZT  pins  are  used.  The  PZT  pins  are 
placed  180  degree  apart  on  a  diameter  of  the  target  plate  holder  and  are  set  flush  with  the 
impact  surface.  The  non-impact  side  of  the  pins  is  connected  directly  to  the  diagnostic 
system  through  a  50Q  coaxial  cables.  A  separate  PZT  pin  is  used  as  the  trigger  pin  for 
activating  the  VISAR  system  as  well  as  the  data  recording  instrumentation. 
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6. 


Details  of  Sample  Preparation 


a.  Ceramic  Corbit-98 

In  order  to  prepare  a  Corbit  target,  first  the  supplied  square  plate  (50  mm 
by  50  mm)  as  reeeived  from  the  manufaeturer  is  eut  to  fit  in  the  round  hole  of  the  target 
holding  plate  using  a  diamond  hole  saw.  Due  to  the  faet  that  eeramie  is  both  a  super  hard 
and  brittle  material,  the  eutting  proeess  is  a  delieate  and  time-eonsuming  proeess.  In 
order  to  ensure  that  the  brittle  samples  do  not  eraek  towards  the  end  when  the  saw  is 
almost  through,  a  saerifieial  pieee  of  eeramie  is  glued  to  the  bottom.  To  prevent  the 
sample  and  the  saw  bit  from  being  overheated  due  to  the  long  eutting  proeess,  the 
eeramie  sample  must  be  immersed  in  a  basin  of  eoolant. 

After  the  eeramie  sample  has  been  eut,  it  is  then  lapped  on  a  rotating 
wheel  affixed  with  a  miero  eloth  whieh  is  eoated  with  9  pm  diamond  suspension  solution, 
followed  by  a  6pm  diamond  solution  to  aehieve  the  required  surfaee  finish.  The  eeramie 
sample  is  then  polished  on  a  diamond  impregnated  lapping  plate  aequired  from  UHL 
Teehnologies,  and  is  designated  DIABLAP.  The  sample  is  typieally  lapped  to  a  flatness 
of  10  pm.  An  example  of  rotating  wheel  used  for  lapping  of  the  target  sample  is  shown  in 
Figure  23. 
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Figure  23.  Buhler  rotating  wheel  used  for  lapping  of  samples. 


Depending  on  the  type  of  experiment  conducted,  a  thin  sheet  of  reflective 
foil,  such  as  a  stainless  steel  shim  of  less  than  5  |um  in  thickness  may  be  glued  to  the  back 
of  the  target  sample  for  the  purpose  of  improving  the  reflection  of  the  VISAR  laser.  The 
target  sample  is  then  glued  into  the  target  holding  plate  and  made  ready  for  the  impact 
experimentation. 

b.  Ceramet 

Due  to  the  high  tolerance  of  the  sample  provided,  no  lapping  was  required 
to  improve  the  flatness  of  the  target.  Although  no  lapping  was  required,  there  was  still  a 
need  to  roughen  the  surface  measured  by  the  VISAR  due  to  the  near  mirror  finish  of  the 
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sample.  This  roughening  was  aeeomplished  by  lapping  the  VISAR  side  of  the  target 
sample  with  diamond  paste.  Typieally  5  |um  diamond  paste  was  used  for  this  roughening 
proeess. 


7,  Design  of  Experiments 
a.  Edge  Effects 

In  order  to  ensure  that  fundamental  assumptions  of  a  uniaxial  strain 
eonditions  holds  true  for  the  experiments  eonducted,  estimation  must  be  made  for 
releases  from  the  edge  of  a  shoek-compressed  sample.  This  estimation  will  help  to  verify 
if  the  inwards  radial  propagation  of  release  waves  generated  at  the  outer  edge  of  the 
sample  ‘pollute’  the  region  of  1-D  strain.  Due  to  the  fact  that  information  about  the 
sound  speed  in  the  compressed  state  as  a  function  of  stress  is  usually  not  known  for  the 
target  material,  this  it  can  only  be  an  estimate.  Table  2  shows  the  equations  used  for  the 
estimation  of  the  time  of  arrival  of  the  edge  release  to  the  center  of  the  sample  where  Ci 
and  Ct  are  the  sound  speeds  of  the  impactor  and  target,  respectively. 
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Table  2.  Equation  used  for  estimating  edge  release. 


t-—  +  — 

^  ~  Cl  ^  Cl  "v 

where  Cj  >  C^, 


where  Q  < 


b.  Spall  Location 

For  symmetrical  impact,  spall  measurement  experiments  are  usually 
conducted  with  impactor  having  half  the  thickness  of  the  target  sample  to  result  in  the 
spall  occurring  in  the  middle  of  the  target.  In  a  situation  where  asymmetrical  impact  is 
done,  or  the  thickness  of  the  sample  cannot  be  changed,  checks  have  to  be  made  to  ensure 
that  spall  occurs  somewhere  in  the  middle  of  the  target  so  that  proper  measurement  of 
spall  strength  can  still  be  made.  This  estimation  can  be  made  using  simple  x-t  diagrams 
as  shown  in  Figure  24. 
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Figure  24.  x-t  diagram  illustrating  interaction  of  release  wave  in  unsymmetrical  impact 

experiment 


In  order  to  simplify  the  process  of  checking  this  for  every  experiment 
conducted,  a  Matlab  code  was  formulated  (see  Appendix  A). 
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IV.  EXPERIMENTAL  RESULTS 


In  this  section,  results  of  all  the  experiments  conducted  for  the  purpose  of  this 
thesis  will  be  presented.  Data  obtained  along  with  the  relevant  analyses  will  be  presented 
here. 

A,  CERAMIC  CORBIT-98  SHOTS 

As  mentioned  in  Denzel  (2010),  Shot  10  3,  which  was  the  only  shot  performed  in 
NFS  for  the  purpose  of  measuring  spall  strength  of  Corbit-98  was  not  successful.  This  is 
suspected  to  be  due  to  “micro-jetting”  of  the  sample  surface  that  destroyed  the  reflectivity 
of  the  foil  and  thus  prevented  any  spall  signature  from  being  recorded.  One  of  the 
objectives  of  this  thesis  was  the  continuation  of  the  research  on  Corbit-98  and  this  section 
documents  the  shots  performed  to  determine  the  spall  strength. 

1.  Shot  1013:  Ceramic  Corbit-98  Spall  Shot 

The  first  planar  impact  experiment  conducted  for  the  purpose  of  this  thesis  was 
designed  to  measure  the  spall  strength  of  Corbit-98  in  a  way  that  was  similar  to  shot 
10  3.  A  thin  layer  of  steel  foil  (12.7  pm  thick)  was  glued  to  the  back  (VISAR  side)  of 
the  ceramic  sample  for  the  purpose  of  reflecting  the  VISAR  laser.  Planned  velocity  was 
0.25  km/s  and  based  on  the  gas  breech  performance  curve.  The  necessary  pressure  to 
achieve  this  velocity  was  approximately  590  psi.  The  actual  velocity  measurement 
yielded  0.2548  +  0.0028  km/s  with  an  indicated  tilt  of  2.04  mrad. 

The  ceramic  Target  Sample  was  measured  to  be  6.179  mm  thick  in  the  center  with 
approximately  a  tilt  of  30  pm  between  surfaces  of  the  target  holding  plate  and  the  sample 
itself.  A  single  crystal  (z-cut)  sapphire  was  used  as  an  impactor.  The  impactor  was 
backed  with  a  piece  of  4.96  mm  thick  closed  cell  foam.  The  inclusion  of  the  foam  was 
necessary  to  provide  a  very  low  shock  impedance  boundary  at  the  back  of  the  impactor  so 
that  a  deep  release  wave  can  be  created  for  the  purpose  of  spalling  the  sample. 
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Figure  25.  Wave  profile  of  ceramic  Corbit-98  spall  shotl0_13 


From  Figure  25,  no  clear  spall  signature  could  be  observed  and  it  is  speculated  to 
be  again  due  to  the  foil  ejecting  off  the  ceramic  target  sample  as  mentioned  by  Denzel 
(2010).  Based  on  theoretical  calculation,  the  approximate  particle  velocity  was  expected 
to  be  0.131  km/s  but  this  was  not  observed  from  the  wave  profile  generated  by  the 
VISAR  system.  Nonetheless,  shock  transit  time  through  the  sample  could  still  be 
determined  and  a  theoretical  particle  velocity  could  still  be  derived  from  the  projectile 
velocity  measured  by  the  velocity  pins. 
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2,  Shot  10  21:  Ceramic  Corbit-98  Spall  Shot 

Due  to  the  possibility  of  the  steel  foil  being  ejeeted  off  the  baek  of  the  target 
sample  in  the  previous  eeramie  spall  shot,  another  spall  shot  was  setup  without  the  use  of 
the  foil.  In  order  to  ensure  that  suffieient  light  was  refleeted  from  the  baek  of  the  target 
sample  to  provide  VISAR  signal,  the  sample  target  was  further  lapped  with  a  6  pm 
diamond  suspension  solution  on  a  Buehler  meohanieally  driven  polisher.  Even  after 
being  polished,  the  VISAR  laser  had  to  be  set  to  twiee  the  usual  wattage  before 
reasonable  signal  amplitude  eould  be  obtained. 

The  planned  projeetile  veloeity  was  0.2  km/s  and  the  required  breeeh  pressure  to 
drive  the  projeetile  is  ealeulated  to  be  388  psi  based  on  the  weight  of  the  projeetile.  The 
eeramie  sample  was  measured  to  be  6.169  mm  thiek  in  the  eenter  with  an  approximate  tilt 
of  30  pm  between  surfaees  of  the  target  holding  plate  and  the  sample  itself.  A  single 
erystal  (z-eut)  sapphire  was  used  as  the  impaetor.  An  air  poeket  was  ereated  at  the  baek 
of  the  impaetor  instead  of  the  foam  used  in  the  previous  shot  for  the  purpose  of  ereating  a 
very  low  shook  impedanoe  boundary. 
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Figure  26.  Wave  profile  of  eeramie  Corbit-98  spall  shotl0_21 


From  Figure  26,  no  elear  spall  signature  could  be  observed  and  peak  particle 
velocity  observed  was  a  low  of  0.03  to  0.04  km/s.  As  this  is  a  symmetrical  impact 
experiment,  the  expected  Ufs  should  be  0.2  km/s.  This  unusually  low  particle  velocity 
and  the  lack  of  spall  signature  indicated  that  the  experiment  has  failed.  The  possible 
reason  of  this  unsuccessful  shot  could  be  due  to  micro  jetting,  due  to  the  shockwave 
created  by  the  planar  impact  hitting  the  imperfect  ceramic  surface  caused  micro  jets  to 
form  and  thus  resulted  in  unusual  data  recorded  by  the  VISAR  system. 

The  flyer  speed  measured  by  the  velocity  pins  was  0.204  +  0.004  km/s  with  a  tilt 
of  4.067  +  2.256  mrad.  Due  to  the  slow  flyer  speed,  the  timing  setup  for  the  oscilloscope 
was  too  short,  and  so  the  time  of  impact  could  not  be  recorded.  ‘Both9.txt’  (see 
Appendix  B)  was  the  file  used  for  this  experiment  and  was  later  modified  to  ‘Bothl  l.txt’ 
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(see  Appendix  C)  in  order  to  be  able  to  reeord  the  impaet  time  for  the  subsequent 
experiment  eondueted  at  this  low  flyer  speed.  Without  the  time  of  impaet,  shoek  speed 
and  partiele  veloeity  eannot  be  determined  for  this  experiment. 

3.  Shot  10  23:  Ceramic  Corbit-98  Spall  Shot 

With  the  unsueeessful  trial  of  experimenting  without  the  use  of  foil,  Shot  10_23 
was  setup  again  using  a  foil  at  the  VISAR  side  of  the  target  sample.  The  foil  used  was  0.3 
pm  made  of  stainless  steel.  In  order  to  ensure  enough  light  was  refleeted  baek  to  the 
VISAR  probe,  the  foil  was  roughened  up  with  a  10  pm  diamond  paste. 

The  planned  veloeity  was  0.2  km/s  and  the  required  breeeh  pressure  to  drive  the 
projeetile  is  ealeulated  to  be  389  psi  based  on  the  weight  of  the  projeetile  measured.  The 
eeramie  sample  was  measured  to  be  6.165  mm  thiek  in  the  eenter  with  an  approximate  tilt 
of  30  pm  between  surfaees  of  the  target  holding  plate  and  the  sample  itself.  A  single 
erystal  (z-eut)  sapphire  was  used  as  the  impaetor.  An  air  poeket  was  ereated  at  the  baek 
of  the  impaetor  instead  of  the  foam  used  in  the  previous  shot  for  the  purpose  of  ereating  a 
very  low  shook  impedanoe  boundary. 

The  impaet  time  of  the  projeetile  was  suooessfully  measured  using  the  flush  pins. 
The  flyer  speed  measured  by  the  veloeity  pins  was  0.2104  ±0.001  km/s  with  a  tilt  of 
1.745  ±  0.6854  mrad.  On  this  last  and  final  experiment  eondueted  to  investigate  the  spall 
strength  of  Corbit-98,  the  VISAR  generated  wave  profde  again  failed  to  reeord  any  spall 
signature  as  shown  in  Figure  27. 
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Figure  27.  Wave  profile  of  ceramic  Corbit-98  spall  shotl0_23 


B.  CERAMET  TUNGSTEN  CARBIDE  SHOTS 

An  investigation  into  a  new  super  hard  material  is  also  an  objective  of  this 
research.  The  new  material  in  focus  for  this  section  of  the  thesis  is  a  metal-ceramic 
composite  material  generally  referred  to  as  ceramet  (Tungsten  Carbide  in  a  Cobalt 
matrix).  Flere  we  tested  two  different  kinds  of  this  material,  called  GC-330  and  GC-915. 
Three  shots  were  performed  on  the  more  dens  GC-915  material  and  of  the  3  shots,  2  were 
designed  to  be  a  free  surface  shot  and  1  was  designed  to  be  a  window  shot.  As  for  GC- 
330,  only  1  free  surface  shot  was  performed.  Appendix  D  includes  the  specifications 
sheets  for  these  materials. 
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1. 


Shot  1014:  Ceramet  Tungsten  Carbide  Spall  Shot 


The  first  Ceramet  shot  condueted  was  designed  to  be  a  free  surface  shot  for  the 
purpose  of  determining  the  spall  strength  as  well  as  obtaining  information  on  the  shock 
Hugoniot  of  this  material.  This  shot  was  setup  as  a  symmetrical  planar  impact  with  the 
impactor  having  a  thickness  of  4.00  mm  and  the  target  having  a  7.996  mm  thickness. 
The  diameters  were  identical,  measuring:  50  mm.  With  symmetrical  impact  as  well  as 
the  selected  thickness  of  impactor  and  target  sample,  the  spall  is  predicted  to  occur  close 
to  the  middle  of  the  target  sample. 

Diagnostics  used  for  these  experiments  include  velocity  pins,  PZT  pins  and  the 
VISAR  systems.  Due  to  the  mirror  like  reflectivity  of  the  Ceramet  target  sample,  there 
was  a  need  to  roughen  the  surface  for  a  more  diffuse  reflection  of  the  laser  beam  from  the 
VISAR.  This  is  necessary  to  ensure  that  the  VISAR  probe  can  receive  the  reflected  laser 
for  the  purpose  of  recording  the  free  surface  motion  of  the  target  sample.  The  roughening 
of  the  surface  was  done  with  the  use  of  5um  diamond  paste  as  shown  in  Figure  28. 


Figure  28.  Diamond  paste  used  for  roughening  of  sample  surface 


Both  surfaces  of  the  target  and  impactor  were  measured  to  have  a  tilt  of  <5  pm 
between  surfaces.  Just  as  for  the  ceramic  spall  shot,  an  air  pocket  was  created  right 


43 


behind  the  impactor  for  the  purpose  of  providing  a  very  low  shock  impedance  boundary 
at  the  back  of  the  impactor  so  that  a  deep  release  wave  is  reflected.  The  mass  of  the 
projectile  with  the  impactor  and  with  the  0-rings  mounted  was  552.9  g.  Using  the  gas 
breech  performance  curve,  the  calculated  pressure  to  achieve  a  desired  projectile  velocity 
of  0.2  km/s  was  446  psi. 

From  the  data  recorded  by  the  velocity  pins,  the  tilt  was  derived  to  be  1.95  + 
0.182  mrad  and  a  projectile  velocity  of  0.209  +  0.0002  km/s.  Shock  transit  time  based 
on  the  impact  time  of  the  flush  PZT  pins  and  the  first  arrival  of  shock  recorded  by  the 
VISAR  system  was  1.091  p,s.  Due  to  offsets  between  the  target  and  the  flush  PZT  pins, 
the  height  differences  must  be  taken  into  account  and  the  corrected  shock  transit  time  was 
1.212  |um.  The  shock  speed  was  calculated  to  be  6.599  km/s.  As  this  is  a  symmetrical 
impact,  the  particle  velocity  Up  can  be  taken  to  be  0.5  *  (Projectile  Velocity)  or 
0.105km/s. 
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Finally,  the  approximate  spall  strength  was  caleulated  to  be  2.11  GPa  using  the 
pull  baek  spall  signature  shown  in  Figure  29.  As  mentioned  earlier  in  the  Chapter  III,  we 
eannot  preeisely  determine  the  spall  strength  from  use  of  bulk  sound  speed  as  the  wave 
speed  in  this  ealeulation  for  brittle  materials.  However,  for  simplifieation,  we  knowingly 
use  it  to  aequire  an  approximate  value  of  spall  strength. 

2.  Shot  1018:  Ceramet  Tungsten  Carbide  Spall  Shot 

A  second  spall  experiment  was  also  performed  to  verify  the  results  of  the  previous 
spall  shot  and  at  the  same  time,  obtain  another  data  point  on  the  shock  Hugoniot. 
Experimental  setup  was  exactly  as  per  described  above  for  the  previous  ceramet  spall 
shot  10_14  with  the  exception  of  a  higher  projectile  velocity.  The  mass  of  the  projectile 
with  the  impactor  and  with  the  0-rings  mounted  was  552.9  g.  From  the  gas  breech 
performance  curve,  the  calculated  pressure  to  achieve  the  desired  projectile  velocity  of 
0.3  km/s  was  1010  psi. 

Wave  profile  generated  by  the  VISAR  system  was  of  very  high  quality.  The  spall 
signature  could  be  clearly  observed  as  shown  in  Figure  30. 


45 


Spall 

Signature 


Figure  30.  Shot  10_18  wave  profile  featuring  spall  signature 


From  the  data  recorded  by  the  velocity  pins,  the  tilt  was  derived  to  be  1 .326  + 
0.37  mrad  and  a  projectile  velocity  of  0.300  +  0.001  km/s.  Shock  transit  time  based  on 
the  impact  time  from  flush  PZT  pins  and  the  first  arrival  of  shock  recorded  by  the  VISAR 
system  was  1.1565  ps.  Due  to  imperfections  from  mounting  the  target  and  PZT  pins, 
height  differences  must  be  taken  into  account  and  the  corrected  shock  transit  time  was 
1.2039  ps.  From  the  target  thickness  and  the  transit  time,  the  shock  was  calculated  to  be 
travelling  at  6.639  km/s.  As  this  is  a  symmetrical  impact,  the  particle  velocity  Up  can  be 
taken  to  be  0.5  (Projectile  Velocity)  or  0.1504  km/s.  Finally,  the  approximate  spall 

strength  was  calculated  to  be  1.959  GPa  using  the  equation  .  As 


mentioned  earlier  in  the  Chapter  III,  we  caimot  precisely  determine  the  spall  strength 
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from  use  of  bulk  sound  speed  as  the  wave  speed  in  this  ealeulation  for  brittle  materials. 
However,  for  simplifieation,  we  knowingly  use  it  to  acquire  an  approximate  value  of 
spall  strength. 


3.  Shot  1019:  Ceramet  Tungsten  Carbide  Window  Shot 

As  opposed  to  the  previous  two  shots  conducted  for  the  ceramet,  Shot  1018  was 
designed  to  be  a  window  shot.  With  the  use  of  a  window,  an  almost  in-situ  measurement 
of  the  particle  velocity  history  could  be  obtained.  A  sapphire  window  was  used  in  this 
experiment  and  as  before,  the  surface  observed  by  the  VISAR  system  was  roughened 
using  6|um  diamond  paste  before  gluing  the  sapphire  window  onto  the  back  of  the  target 
sample. 

In  order  to  measure  the  HEL  of  this  ceramet,  a  higher  projectile  velocity  is  desired 
to  increase  the  impact  pressure.  The  experiment  was  setup  to  throw  the  projectile  at  0.35 
km/s.  In  order  to  reach  this  velocity  without  using  exceedingly  high  breech  pressure,  the 
generic  projectile  normally  used  was  bored  out  internally  to  lower  the  mass  to  447.6  g. 
The  required  breech  pressure  was  estimated  to  be  1,100  psi  with  the  use  of  the  gas  gun 
performance  curve  to  achieve  this  required  projectile  velocity. 

This  shot  was  setup  as  a  symmetrical  planar  impact  with  the  impactor  measured  at 
4.00  mm  thickness  at  the  center  and  the  target  measuring  7.993  mm  thickness  at  the 
center — the  diameters  were  the  same  as  before.  With  symmetrical  impact,  the  Up  can  be 
taken  to  be  0.5  (Projectile  Velocity).  Both  surfaces  of  the  target  and  impactor  were 

measure  to  have  a  tilt  of  <5  pm  between  surfaces. 
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Figure  3 1 .  Shot  1019  wave  profile  of  a  window  shot 


All  diagnostics  used  correctly  recorded  excellent  data.  As  shown  in  Figure  31,  the 
wave  profile  generated  by  the  VISAR  system  was  of  very  high  quality  and  a  subtle 
precursor  was  observed  at  the  toe  of  the  wave  profile.  Based  upon  the  shorting  time  of 
the  velocity  pins,  the  projectile  was  calculated  to  be  travelling  at  0.35  +  0.001  km/s  and 
had  a  tilt  of  2.987  +  0.604  mrad.  Taking  into  account  the  imperfections  in  mounting  of 
the  target  sample  and  the  PZT  fiush  pins,  the  shock  transit  time  was  corrected  to  1.080  qs 
and  the  resulting  shock  speed  Us  was  calculated  to  be  7.402  km/s.  The  particle  velocity 
Up  was  taken  to  be  0.5  (Projectile  Velocity)  or  0.175  km/s. 
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4,  Shot  10  20:  GC-330  Ceramet  Tungsten  Carbide  Spall  Shot 

This  final  experiment  condueted  on  a  eeramet  was  done  to  eompare  the  shoek 
Hugoniot  of  a  less  dense  ceramet  (GC-330).  The  experiment  was  designed  to  be  a  free 
surface,  symmetrical  planar  impact.  The  experimental  setup  was  close  to  shot  1018 
with  the  exception  of  a  lower  breech  pressure  (990  psi)  used  for  this  experiment. 
Although  the  projectile  is  fired  at  the  same  velocity  of  0.3  km/s,  a  lower  pressure  was 
required  due  to  a  less  dense  impactor,  which  brings  down  the  mass  of  the  projectile. 


Figure  32.  Shot  10_20  wave  profile  of  a  window  shot 


As  shown  in  Figure  32,  the  wave  profile  generated  by  the  VISAR  system  is  of 
very  high  quality  and  spall  signature  can  be  clearly  observed.  From  the  velocity  pins,  the 
projectile  was  calculated  to  have  a  velocity  of  0.318  +0.001  km/s  and  a  tilt  of  1.241 
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±  0.366  mrad.  Shock  transit  time  was  measured  to  be  1.218  |us  but  eorreeted  to  1.287  |us 
after  taking  into  aceount  the  imperfeetion  of  the  mounting  of  the  target  and  PZT  pins. 
The  shoek  speed  Us  was  ealeulated  to  be  6.216  km/s  based  on  the  transit  time  and 
thiekness  of  the  target.  The  partiele  veloeity  Up  ean  be  taken  to  be  0.5  *t/^(Projeetile 
Veloeity)  or  0.159  km/s  due  to  symmetrieal  impaet.  Finally,  the  approximate  spall 
strength  was  ealeulated  to  be  2.284  GPa. 
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V.  DISCUSSION  OF  RESULTS 


A,  UNCERTAINTY  ANALYSIS 

In  order  to  understand  how  aeeurate  the  measurements  deseribed  above  are, 
there’s  a  need  to  look  at  the  potential  sourees  of  uneertainty  and  make  an  estimate  using 
standard  teehniques.  In  this  seetion,  the  prineipal  sourees  of  error  are  measured  in 
distances  and  times. 

Due  to  the  design  of  the  gas  gun  (throwing  projectile  at  very  low  speed),  a  much 
higher  shock  speed  exists  in  the  target  sample  as  compared  to  the  projectile  velocity,  and 
this  results  a  relatively  short  shock  transit  time.  This  results  in  a  large  uncertainty  when  it 
comes  to  establishing  shock  transit  time,  because  small  distance  offsets  in  flush  pins  lead 
to  relatively  large  time  offsets.  With  the  assumption  of  a  Gaussian  distribution  of 
uncertainty,  fractional  uncertainty  in  shock  velocity  can  be  calculated  based  on  the 
equation  below: 


Although  the  electronic  measuring  equipment  in  the  lab  allows  for  measurement 
of  up  to  1-2  qm,  the  limiting  factor  lies  in  the  flatness  of  the  sample  target  and  how 
parallel  the  front  and  back  surfaces  are  to  each  other.  Therefore,  Ax  is  based  on  the 
flatness  tolerance  measured  where  x  is  the  thickness  of  the  sample.  At  was  set  at  15  ns 
which  is  the  resolution  of  time  the  diagnostics  can  measure  to  and  t  is  the  transit  time 
measured.  Table  3  summarizes  the  values  used  for  the  calculations  of  uncertainty  in 
measured  shock  velocity. 
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Table  3.  Values  used  for  ealeulating  uneertainty 


Ceramic 

Ceramet 

Shotl0_13 

Shotl0_21 

Shotl0_23 

Shotl0_14 

Shotl0_18 

Shotl0_19 

Shotl0_20 

Ax  (mm) 

0.010 

0.015 

0.0013 

0.003 

0.003 

0.003 

0.003 

X  (mm) 

6.179 

6.169 

6.165 

7.996 

7.993 

7.993 

8.003 

Ax 

—  (%) 

X 

0.16 

0.24 

0.02 

0.04 

0.04 

0.04 

0.04 

At  (ns) 

15 

15 

15 

15 

15 

15 

15 

t(|is) 

0.571 

- 

0.573 

1.211 

1.204 

1.08 

1.287 

o  1 

2.63 

- 

2.62 

1.24 

1.25 

1.39 

1.17 

At 

t 

2.63 

- 

2.62 

1.24 

1.25 

1.39 

1.17 

Another  potential  souree  of  uneertainty  lies  in  the  eorrect  ealeulation  of  shook 
stress.  With  the  assumption  of  very  good  measurement  of  initial  density  of  sample  target 
used,  uneertainty  of  shook  stress  oan  be  easily  oaloulated  with  Equation  21: 


AP  ( AUX  ( AUX 


(21) 


The  uneertainty  of  in  impaot  velooity  was  determined  through  the  use  of  a  least  square  fit 
on  the  shorting  time  of  the  velooity  pins.  Using  the  uneertainty  of  the  shook  speed  derived 
in  earlier,  the  oaloulated  uncertainties  for  the  respective  shock  stresses  are  tabulated  in 
Table  4. 
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Table  4.  Values  used  for  ealculating  shoek  stress  uneertainty 


Ceramic 

Ceramet 

Shotl0_13 

Shotl0_21 

Shotl0_23 

Shotl0_14 

Shotl0_18 

Shotl0_19 

Shotl0_20 

AU, 

- ^(%) 

U 

2.63 

- 

2.62 

1.24 

1.25 

1.39 

1.17 

AU^ 

- -  (%) 

u 

1.10 

1.96 

0.48 

0.96 

0.33 

0.35 

0.31 

AP 

—  (%) 

2.63 

- 

2.62 

1.24 

1.25 

1.39 

1.17 

B,  CEARMIC  CORBIT-98 

1,  Spall  Strength 

As  a  eontinuation  of  previous  NPS  researeh  performed  by  Denzel  (  2010),  the 
main  foeus  on  eeramie  Corbit-98  was  to  establish  and  verify  the  spall  strength  as  well  as 
adding  more  data  points  on  the  shoek  Hugoniot  in  the  Us  vs  Up  plane.  Unfortunately,  all 
free  surfaee  planar  impaet  experiment  eondueted  on  the  NPS  gas  gun  untill  the  present 
date  has  failed  to  reeord  any  elear  spall  signature.  Mentioned  earlier  in  Chapter  IV, 
jetting  from  the  porous  surfaee  of  the  eeramie  sample  is  the  most  likely  eause  of  this. 

2.  Hugoniot 


Tables.  Hugoniot  data  of  Corbit-98 


Shot 

Material 

Type 

Ud  (km/s) 

Us  (km/s) 

Up  (km/s) 

P  (GPa) 

1013 

Corbit  98 

FS 

0.254 

10.81 

0.1325 

5.515 

1021 

Corbit  98 

FS 

0.204 

- 

- 

- 

10_23 

Corbit  98 

FS 

0.210 

10.756 

0.1328 

5.499 
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Values  for  the  elastie  Hugoniot  for  the  ceramic  shots  conducted  for  the  purpose  of 
this  thesis  are  summarized  in  Table  5.  Due  to  the  free  surface  micro-jetting  effects 
described  above,  Up  as  measured  from  the  wave  profde  cannot  be  accurately  used  for  the 
calculation  of  the  pressure  state  attained.  Therefore,  Up  was  calculated  based  on 
impedance  matching  method  using  the  measured  Ud  and  Us-  In  order  to  simply  the 
process,  a  simple  Matlab®  code  (see  Appendix  E)  was  generated  for  the  purpose  of  this 
calculation. 


Figure  33.  Compiled  Hugoniot  data  for  ceramic  (After  Denzel,  2010) 

The  elastic  Hugoniot  for  Ceramic  is  particularly  difficult  to  measure  due  to  the 

fact  that  a  slow  projectile  velocity  with  high  shock  velocity  in  target  sample  can  result  in 

large  time  uncertainties  with  just  a  small  distance  uncertainties.  From  Figure  33,  it  can 

be  observed  that  the  current  data  are  in  reasonable  agreement  with  the  Hugoniot 
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previously  measured  by  Denzel  (2010),  changing  only  slightly  from  the  previous 
=  10.53  +  0.936f/p  to  the  current  derived  Hugoniot  of  =  10.57  +  0.95917/^ .  With 

the  inclusion  of  the  two  more  data  points  from  the  experiment  of  this  thesis,  the  new 
elastic  Hugoniot  found  is  a  good  starting  point  for  computer  simulation  of  Corbit-98. 

C.  CERAMET  TUNGSTEN  CARBIDE 
1.  Hugoniot 


Table  6.  Hugoniot  data  obtained  for  Ceramet  Tungsten  Carbide 


Shot 

Material 

Type 

Ud 

(km/s) 

Us  (km/s) 

Up  (km/s) 

P  (GPa) 

1014 

GC-915 

FS 

0.209 

6.599 

0.105 

9.276 

1018 

GC-915 

FS 

0.301 

6.639 

0.150 

14.006 

1019 

GC-915 

SW 

0.350 

7.402 

0.175 

18.154 

1020 

GC-330 

FS 

0.318 

6.216 

0.159 

12.351 

Table  6  shown  above  summarizes  all  the  Hugoniot  data  collected  for  the  Ceramet 
Tungsten  Carbide  experimented  in  support  of  this  thesis.  Figures  tabulated  are  already 
corrected  for  the  imperfections  of  the  mounting  of  the  sample  and  PZT  pins.  Using  the 
data  in  Table  3,  the  Us-U?  plot  is  as  shown  in  Figure  34  and  the  linear  fit  is  based  solely 
on  GC-915  data  only.  The  GC-330  data  point  is  shown  for  completeness  only. 
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Figure  34.  Flugoniot  plot  of  NFS  ceramet  data 

2,  Dynamic  Strength 

One  of  the  main  objectives  of  this  research  was  to  determine  the  compressive 
dynamic  yield  strength  for  Ceramet  Tungsten  Carbide.  In  order  to  do  this,  there  is  a  need 
to  shock  the  sample  to  a  pressure  state  above  the  yield  point.  Due  to  velocity  limitations 
of  the  single  stage  gas  gun  used  for  the  purpose  of  this  research,  the  maximum  stress  state 
achievable  is  not  significantly  higher  than  the  HEL.  Figure  35  shows  the  plot  of  the 
compression  shock  profile  of  the  ceramet  experiment  conducted  for  this  research; 


56 


Figure  35 .  Compression  shock  profiles  for  tungsten  carbide 


From  Figure  35,  it  can  be  observed  that  the  compression  shock  profile  of  NFS 
does  not  display  a  distinct  transition  from  the  elastic  precursor  rise  to  the  transition  ramp 
region.  This  phenomenon  was  also  observed  by  Grady  (1995). 
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TIME  (ms) 

Figure  36.  Compression  shook  profiles  for  tungsten  carbide  (From  Grady,  1995) 

From  Figure  36,  it  can  be  observed  that  WC-3  compression  shock  profile  did  not 
provide  a  distinct  transition  level  as  WC-1  and  WC-2.  It  is  noteworthy  that  both  WC-3 
and  NFS  compression  shock  profiles  are  at  a  low  peak  Up  of  0.2  to  0.3  km/s.  By  firing  at 
higher  impact  speeds,  larger  shock  amplitudes  will  be  reached  that  should  provide  a  more 
distinct  yielding  transition  as  those  observed  in  Grady’s  WC-2  and  WC-1  profiles  in 
Figure  35. 

There  is  a  consistent  transition  point  at  approximately  0.01  km/s  that  indicates  the 
start  of  the  yielding  process  (i.e.,  HEL  stress).  Based  on  a  shock  impedance  calculation 
this  leads  to  a  HEL  of  0.935  GPa.  However,  we  believe  this  is  just  the  beginning  point  for 
a  complex  yielding  process  that  is  spread  out  over  a  broad  range  of  particle  velocity  and 
stress.  In  our  wave  profiles,  we  see  one  and  perhaps  two  changes  in  slope  (slightly  above 
Up  of  0.1  km/s  and  a  more  distinct  change  slightly  below  0.15  km/s).  These  may  suggest 
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that  the  HEL  point  at  0.01  km/s  indicates  the  onset  of  a  partial  yielding  process  and  the 
0.15  km/s  point  shows  a  full  yielding  process.  The  intricacies  of  this  complex  yielding 
are  far  beyond  the  scope  of  this  thesis  and  more  data  will  be  needed  to  obtain  a  better 
understanding  of  this  process.  We  have  been  collaborating  with  Prof  Zok  at  UCSB  and 
one  of  his  PhD  students,  Mr.  Brett  Compton,  will  be  investigating  this  complex  yielding 
phenomena. 


3,  Spall  Strength 


Figure  37.  Wave  profdes  of  NPS  tungsten  carbide  spall  shots 
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Clear  spall  signatures  can  be  easily  observed  for  all  the  wave  profiles  obtained  as 
shown  in  Figure  37.  These  clear  signatures  indicated  that  the  material  has  completely 
spalled.  The  spall  strength  of  tungsten  carbide  can  then  be  calculated  from  the  spall 
pullback  signal  identified  in  the  compression  shock  profile. 

Similar  to  ceramic  materials,  it  is  not  rigorous  to  use  bulk  sound  speed  to 
calculate  the  spall  strength  as  would  be  done  for  metal  that  is  of  a  more  ductile  nature. 
Therefore  the  spall  strength  calculated  here  is  only  an  approximation  done  for 
comparison,  and  a  small  error  will  exist  in  our  calculated  values.  A  more  rigorous 
analysis  of  these  data  is  also  beyond  the  scope  of  this  thesis. 


Table  7.  Spall  properties  of  tungsten  carbide 


Shot 

Material 

^max  (km/s) 

^min  (km/s) 

(^spaii  (GPa) 

10-14 

GC-915 

0.2003 

0.139 

2.11 

10-18 

GC-915 

0.289 

0.232 

1.96 

10-20 

GC-330 

0.31 

0.234 

2.28 

Table  7  summaries  spall  strength  calculated  from  the  present  data.  The  bulk  sound 
speeds  used  in  the  above  calculation  are  4.69  km/s  and  4.9  km/s  for  GC-915  and  GC-330 
respectively.  A  reasonable  agreement  with  the  spall  strength  of  between  2.62-3.56  GPa 
as  reported  by  Grady  (1995). 
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VI.  CONCLUSIONS 


From  the  experiments  performed  herein,  fundamental  dynamic  loading  properties 
for  both  ceramet  tungsten  carbide  GC-915  and  GC-330  were  found.  A  summary  of  the 
results  are  tabulated  in  Table  8. 


Tables.  Summary  of  results 


Shot 

Material 

Type 

Ud  (km/s) 

Us  (km/s) 

Up  (km/s) 

^ spall  (GPa) 

P  (GPa) 

1013 

Corbit  98 

ES 

0.254 

10.81 

0.1325 

- 

5.515 

10_21 

Corbit  98 

ES 

0.204 

- 

- 

- 

- 

10_23 

Corbit  98 

ES 

0.210 

10.756 

0.1328 

- 

5.499 

1014 

GC-915 

ES 

0.209 

6.599 

0.105 

2.11 

9.276 

1018 

GC-915 

ES 

0.301 

6.639 

0.150 

1.96 

14.006 

1019 

GC-915 

SW 

0.350 

7.402 

0.175 

N.A. 

18.154 

1920 

GC-330 

ES 

0.318 

6.216 

0.159 

2.28 

12.351 

With  3  planar  impact  experiments  conducted  for  GC-915  tungsten  carbide,  there 
is  enough  information  on  the  Hugoniot  in  Us-Up  plane  to  best  fit  a  Hugoniot  relationship 
as  follows: 

f/5  =  10.2  + 5.42 

Since  the  Hugoniot  relationship  for  tungsten  carbide  is  only  derived  based  on  3  points, 
there’s  still  a  fair  amount  of  uncertainty  and  a  few  more  experiments  are  recommended  to 
verify  the  reliability.  Above  all,  the  experiments  conducted  for  the  purpose  of  this  thesis 
were  only  performed  at  relatively  low  impact  velocity  due  to  the  limitation  of  the  NFS 
gas  gun.  Although  these  low  impact  velocities  give  us  a  good  preview  of  the  dynamic 
loading  response  of  the  material  tested,  it  is  necessary  to  conduct  similar  experiments  at  a 
higher  stress  state  to  obtain  fast  rising  plastic  waves  and  be  able  to  better  understand  the 
yielding  process. 

With  regards  to  the  ceramic,  we  were  able  to  obtain  additional  Hugoniot  data  that 


reinforces  the  EOS  obtained  by  Denzel  (2010).  Our  Hugoniot  EOS  changes  only  slightly 
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from  the  previous  =  10.53  +  0.936f/^  to  the  eurrently  derived  Hugoniot  of 
=  10.57  + 0.959 If/^.  Due  to  miero  jetting  from  the  porous  surfaee  of  the  target 

sample,  the  generated  VISAR  wave  profdes  eould  not  be  used  to  determine  the  spall 
strength.  More  investigation  into  the  experimental  method  for  establishing  spall  strength 
is  reeommended  for  future  researeh  but  is  beyond  the  seope  of  this  researeh. 

Corbit-98  eeramie  has  low  spall  strength  of  2.71  kbar  (from  only  1  experiment 
eondueted  by  Denzel  (2010))  illuminates  the  problem  that  alumina  based  eeramies  have 
poor  or  no  multi-hit  eapability  for  armor  applieations.  Espeeially  for  the  eomposite 
armor  eoneept,  the  seeond  wave  spreading  layer  that  is  usually  made  of  low  impedanee 
material,  will  refleet  a  release  wave  baek  into  the  eeramie  layer  eausing  it  to  spall. 
Although  eeramet  tungsten  earbide  has  a  spall  strength  of  almost  10  times  higher  than 
Corbit-98  (2.11  GPa  and  0.27  GPa,  respeetively),  it  is  still  relatively  low,  and  the  ideal 
first  layer  for  the  eomposite  armor  eoneept  would  be  a  material  with  high  HEL  and  spall 
strength  but  sueh  materials  does  not  eurrently  exist.  Promising  results  have  been  shown  in 
researeh  of  amorphous  metals  but  these  materials  are  also  known  to  be  relatively  brittle. 
This  inherent  brittle  response  ean  be  improved  by  having  duetile  dendritie  inelusions 
gives  amorphous  metals  the  ability  to  limit  the  distanees  a  eraek  ean  propagate  and  thus 
inereasing  the  overall  material’s  duetility.  It  is  reeommended  that  future  work  on  planar 
impaet  experiments  be  eondueted  on  amorphous/dendritie  metals  to  better  eharaeterise 
their  dynamie  loading  response.  Eor  now,  small  tiles  of  eeramie  or  even  eeramet  ean  be 
assembled  to  form  the  first  layer  in  bid  to  obtain  a  better  multi-hit  eapability  for  the 
eomposite  armor  eoneept. 

With  the  new  information  aequired  in  thesis,  simulation  ean  now  be  improved  for 
use  in  hydrodynamie  eodes.  Using  Johnson-Holmquist  model  with  better  input,  a  new 
eomputer  simulation  ean  be  done  based  on  Ongs  (2009)  work.  Eurthermore,  new 
information  on  eeramet  tungsten  earbide  is  now  available  to  be  tested  in  the  simulation 
model  of  the  eomposite  armor  system.  These  simulations  ean  then  be  used  to  eompare 
with  existing  and  new  integral  test  results. 
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APPENDIX  A:  X-T  DIAGRAMS  MATLAB  CODE 


%%  Initial  Condition 
clc 

clear  all 

rho2=3.85;  %Ceramic 

rhol=3.985;  %Sapphire 

cl=ll . 19; 

c2=l  0  .  5  ; 

sl  =  l; 

s2=0 . 94; 

Ud=0 .25; 
xl=-4 : 0; 
x2=0 : 6; 

%%  Coefficient  of  quadratic  equation 
a=(rhol*sl-rho2*s2) 
b=-  (2*rhol*sl*Ud+rhol*cl  +  rho2*c2) 
c=(rhol*sl*Ud'^2  +  rhol*cl*Ud) 

%%  Solve  quadratic  equation 

Upxl=  -  (b  +  (b'"2  -  4*a*c)  (1/2)  )  /  (2*a) 

Upx2=  -  (b  -  (b'"2  -  4*a*c)  (1/2)  )  /  (2*a) 

%%  Ask  for  user's  choice  of  Up2 
Up2=input ( ' Upxl  or  Upx2  for  Up2?  ') 

%%  Us2 

Us2=10 . 5+0 . 94*Up2 
%%  First  wave 

t2=x2/Us2  %wave  in  ceramic 

plot  (x2 , t2 ,  ' b ' ) 
hold  all 

Usl=- ( (Ud-Up2) +cl) ; 

tl=xl/Usl  %wave  in  sapphire 

plot  (xl, tl,  ' r ' ) 

%%  Reflected  wave 

t3=x2/ (1.2*Us2)  %reflected  wave  in  ceramic 

plot  (x2,fliplr(t3) +max (t2 ) , ' --b ' ) 
t4=xl/ (1 .2*Usl) 

plot  (xl,fliplr(t4) +max (t 1 ) , ' --r ' ) 

plot  (x2,t2*l. 2+max (fliplr (t4) +max (tl)), '--r') 

Timeofspall=max (tl) +4/ (1.2*Us2) 
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APPENDIX  B:  BOTH9.TXT 


RunDig  setup  file. 

Experiment  description  line 

DEVICE :DP04 104 : BO 201 41 
Digitizer  use  =  VISAR 
Factory 
Select: CHI  ON; 

CHliCoupling  DC; Impedance  50;Pos  +4.0;Scale  0.02; 

Select :CH2  ON; 

CH2:Coupling  DC; Impedance  50;Pos  +4.0;Scale  0.02; 

Select :CH3  ON; 

CH3:Coupling  DC; Impedance  50;Pos  +4.0;Scale  0.02; 

Select :CH4  ON; 

CH4: Coupling  DC; Impedance  meg;Pos  +4.0; Scale  0.02; 

Horiz : RecordLength  100000; Scale  2E-6; delay : time  8.0e-6; 
Trigger : A: Type  Edge; 

Trigger :A:Edge : Source  EXT; Coup  DC; Slope  rise; 

Trigger :A:Mode  Normal;Level  1.00; 

DEVICE:DPO4104 :B020142 
Digitizer  use  =  pins 
Factory 
Select: CHI  ON; 

CHl:Coupling  DC; Impedance  meg;Pos  +3.0;Scale  1.00; 

Select :CH2  ON; 

CH2:Coupling  DC; Impedance  meg;Pos  -4.0; Scale  10.0; 

Select :CH3  ON; 

CH3:Coupling  DC; Impedance  meg;Pos  -4.0; Scale  10.0; 

Select :CH4  ON; 

CH4:Coupling  DC; Impedance  meg;Pos  +0.0;Scale  0.05; 

Horiz : RecordLength  100000;Scale  lE-5 ; delay : time  -lO.Oe-6; 
Trigger : A: Type  Edge; 

Trigger :A:Edge : Source  EXT; Coup  DC; Slope  rise; 

Trigger :A:Mode  Normal;Level  1.00; 
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APPENDIX  C:  BOTHll.TXT 


RunDig  setup  file. 

Experiment  description  line 

DEVICE :DP04 104 : BO 201 41 
Digitizer  use  =  VISAR 
Factory 
Select: CHI  ON; 

CHliCoupling  DC; Impedance  50;Pos  +4.0;Scale  0.02; 

Select :CH2  ON; 

CH2:Coupling  DC; Impedance  50;Pos  +4.0;Scale  0.02; 

Select :CH3  ON; 

CH3:Coupling  DC; Impedance  50;Pos  +4.0;Scale  0.02; 

Select :CH4  ON; 

CH4: Coupling  DC; Impedance  meg;Pos  +4.0; Scale  0.02; 

Horiz : RecordLength  100000; Scale  2E-6; delay : time  8.0e-6; 
Trigger : A: Type  Edge; 

Trigger :A:Edge : Source  EXT; Coup  DC; Slope  rise; 

Trigger :A:Mode  Normal;Level  1.00; 

DEVICE:DPO4104 :B020142 
Digitizer  use  =  pins 
Factory 
Select: CHI  ON; 

CHl:Coupling  DC; Impedance  meg;Pos  +3.0;Scale  1.00; 

Select :CH2  ON; 

CH2: Coupling  DC; Impedance  meg;Pos  -4.0; Scale  10.0; 

Select :CH3  ON; 

CH3:Coupling  DC; Impedance  meg;Pos  -4.0; Scale  10.0; 

Select :CH4  ON; 

CH4:Coupling  DC; Impedance  meg;Pos  +0.0;Scale  0.05; 

Horiz : RecordLength  100000;Scale  2E-5 ; delay : time  -lO.Oe-6; 
Trigger : A: Type  Edge; 

Trigger :A:Edge : Source  EXT; Coup  DC; Slope  rise; 

Trigger :A:Mode  Normal;Level  1.00; 
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APPENDIX  D:  MATLAB  CODE  FOR  SOLVING  Up 


%%  Initial  Condition 
clc 

clear  all 

rho2=3.85;  %Ceramic 
rhol=3.985;  %Sapphire 
cl=ll . 19; 
sl  =  l; 

Ud=0 .2548; 
xl=-4 .018:0; 
x2=0 : 6. 179; 

Us2=14 .369 

Us2=input ( ' Input  the  experiment  Us  found?  ' ) 
%%  Coefficient  of  quadratic  equation 
a=rhol*sl 

b=-  (2*rhol*sl*Ud+rhol*cl  +  rho2*Us2) 
c=(rhol*sl*Ud'^2  +  rhol*cl*Ud) 

%%  Solve  quadratic  equation 

Upxl=  -  (b  +  (b'"2  -  4*a*c)  (1/2)  )  /  (2*a) 

Upx2=  -  (b  -  (b"'2  -  4*a*c)  (1/2)  )  /  (2*a) 

Up2=input ( ' Upxl  or  Upx2  for  Up2?  ') 
P=rho2*Us2*Up2 
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APPENDIX  E:  GRADE  SPECIFICATIONS  FOR  CERAMET 


GEISE9AI  CARBIDE  PRODJCTE  ARE  MADE  lA  THE  USA 


Grade  Specifications 


GENERALCARBIDE. 


Avoiga 


ChamlcBl  Comperfllon  ‘nBimwBM 

Wd^it  Pcimit  HwdiMCB  OvRsIQr  Rupun 


wc 

Cb 

Qdier 

HM 

Sfiuij|lh  (pd) 

Tungsten  Caifaida  Gradas  wU  Oibalt  KndBr 

Ofimlcion  rmbnriooiO 

GG«12f 

86 

12 

922-932 

14X16-1420 

471000 

GC01SF* 

85 

15 

908-918 

1179-1192 

56Cy000 

0.8  iTicion  rsubmioroiO 

GC<n5 

M3 

53 

914-942 

14J88-14^7 

441000 

GGOKP 

90 

10 

914-922 

1439-1431 

520000 

CC^)10CR" 

89 

10 

1 

923-933 

1425-1435 

530000 

GC015* 

85 

15 

893-903 

13B9-14J03 

531000 

GC-015CR* 

84 

15 

1 

904-914 

1324-1186 

550000 

ijOrnfoiMi 

GC-10B 

963 

17 

92J-933 

15.12-1521 

480X100 

GC-106 

94 

6 

919-92J 

14JB6-1497 

510X100 

GC-109 

91 

9 

91JD-918 

1434  - 1436 

525X100 

2j0  macron 

GC-206 

94 

6 

912-922 

1486-1497 

4BOJOOO 

GC-2D9 

91 

9 

902-912 

1433  - 1435 

511000 

GC-211* 

89 

11 

894-904 

1433  - 14.45 

SSOyOOO 

BjOiraaon 

ec-310* 

90 

10 

893-903 

1446-1438 

510X»0 

GC-313* 

87 

13 

88.1-89.1 

14.15-1427 

511000 

GC-31S* 

85 

15 

873-883 

119S-14J09 

520000 

GC-32<r 

80 

20 

853-863 

1346-1334 

470000 

GC-325* 

75 

25 

833-84J 

1103-1123 

440X100 

GC-33(P 

70 

30 

813-829 

1231  - 1282 

420X100 

dOincron 

GC-950* 

85 

15 

834-874 

1195-14J09 

480X100 

GC618* 

82 

18 

852-862 

1337- 1181 

450XX» 

ITOrncran 

QC-siy 

85 

15 

853>863 

1195-1409 

AKWIDO 

ijOni&OmIcnn  Aitad  flfuaun) 

GC-AMM 

94 

6 

904-914 

1489-1500 

480000 

GC-eOBM 

9(2 

8 

B0B-90M 

1436-1428 

490000 

GC-filOM 

90 

10 

888-898 

1446-1438 

500000 

GC^ZM 

SB 

12 

882-882 

1425  1437 

51O000 

SintcrHIP  ProcMS  Giwrantewi 


IISIGvdEnSInxt 
Graembwg,PAlSCOl  USA 


T  ■0aZ45.24M  +i  724A3&3D00 
F  +1  72'U3&j6274 


YHW  prnrniiniihlffrinnfn 
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